Non-Markovian control of qubit thermodynamics by frequent quantum measurements 
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We explore the effects of frequent, impulsive quantum nondemolition measurements of the energy 
of two- level systems (TLS), alias qubits, in contact with a thermal bath. The resulting entropy 
and temperature of both the system and the bath are found to be completely determined by the 
measurement rate, and unrelated to what is expected by standard thermodynamical rules that hold 
for Markovian baths. These anomalies allow for very fast control of heating, cooling and state- 
purification (entropy reduction) of qubits, much sooner than their thermal equilibration time. 



I. INTRODUCTION 



Non-Markovian quantum thermodynamics of two level 
systems (TLS) in contact with a bath has surprising as- 
pects in store. According to standard Markov thermo- 
dynamics, the TLS (alias qubit) thermal equilibration 
process is expected to progress monotonically, accom- 
panied by increase of the entropy, at least on average 
Yet drastic deviations from this trend are 
revealed when considering impulsive disturbances of ther- 
mal equilibrium between TLS and a bath p, @] ■ These 
effects bear certain similarities to the work described in 
We have shown [§| that frequent and brief quan- 
tum non demolition (QND) measurements of the TLS 
energy-states entail unfamiliar anomalies of the entropy 
and temperature of both system and bath, which become 
unrelated to what is known from standard, Markovian 
thermodynamic rules 0, HI]: (i) a transition from heating 
to cooling of the TLS ensemble as we vary the interval be- 
tween consecutive measurements on the time scale of the 
inverse energy separation of the qubit levels; and (ii) cor- 
respondingly, oscillations of the entropy relative to that 
of the equilibrium state. 



II. SYSTEM-BATH ENTANGLEMENT AT 
EQUILIBRIUM 

A. Hamiltonian 



The following Hamiltonian describes the qubit system 
that interacts with the bath. 



H to t — Ha + Hr 



H. 



SB- 



(1) 



Here H tot pertains to the coupled system and bath and 
consists of: 



H s =hio a \e)(e\, 
H B =hy^ j u\a\a\, 



(2) 
(3) 



H S b =SB,S = a x , B = [K\ax + K* x a\) , (4) 



where S and B are the system and bath operators, re- 
spectively, in the system-bath interaction H$b, a \( a \) 
are the annihilation (creation) operators, and k\ is the 
matrix element of the weak coupling to bath mode A. We 
stress that in the interaction Hamiltonian (Hsb) we do 
not invoke the rotating- wave approximation (RWA) (Iol| . 
namely, we do not impose energy conservation be- 
tween the system and the bath, on the time scales 
considered [ill]. 



Here we present an in-depth study of short-time evolu- 
tion of quantum systems coupled to a bath, interrupted 
by frequent measurements. We first discuss in Sec.HTlthc 
initial equilibrium state relevant to our scenario. Sec.lIIII 
then describes the measurement-induced disturbance of 
equilibrium. In Sec. IIVI we present a master equation 
analysis of the post-measurement evolution and a discus- 
sion of the heating and cooling requirements. Cooling 
conditions and entropy evolution of the system are dis- 
cussed in Sec. [V] and |VH respectively. A discussion of 
possible experimental realizations is given in Sec. IVIll 



B. Qubit state mixedness at equilibrium 

At equilibrium, the qubit and the bath are in an en- 
tangled state. To find the mean energy mixedness (im- 
purity) of the qubit (TLS) at a given temperature T, one 
needs the equilibrium density matrix for the total system 
PEq = c:x v{~f3Ht t)/Z , where Z is the partition function 
and 1//3 = k B T. 

Using Heims perturbation theory [l2| one can expand 



2 



PEq 



PEq = 7?e 



-P(H +H S t 



-13H 



z 



[l + eS 1 +e 2 S 2 + 0(e 3 ) + ■■■], (5) 



where 



e = m&x(r) k /fkJ a ), 

is a small dimensionless parameter normalizing the rate 
r/k of the maximally coupled mode to the TLS natural 
frequency, and 



eSx=p[ dxe xf3H "H SB e~ xl3Ho , 
Jo 



(6a) 



The frequency — and temperature — dependent coeffi- 
cients in ([7]) are 



1 



{(cosh (3 h(uj a =F w) — 1) 



±(sinh/5ft(u; Q T^)-^KT^))]- (12) 

From Eq. ([7]) it can be seen that even at zero temper- 
ature purity is incomplete, P Eq (e) < 1, which is due to 
the system-bath entanglement. The difference between 
P Eq in ([9]) and P Eq (e) in ((Jj has a non-monotonic depen- 
dence on p. This can be seen from Fig. [T] where we have 
plotted the relative change of TLS purity with inverse 
temperature. As the purity drop that we wish to cor- 
rect is non-monotonic with temperature, so will be the 
resultant purification. 



e 2 S 2 



(3 I dx I dye 
Jo Jo 



xf3H 



H SB e- (x - v)fiHo HsBe~ yl3Ho - (6b) 



Noting that e~P H ° = Z ps <S> pb, where ps and p E are 
the equilibrium density matrices for the system and the 
bath without interaction, the trace over the bath degrees 
of freedom can be performed. The state of the system is 
diagonal in the a z basis and is given by 

Ps(e) = \(I+P Eq (e)<J z ). 

The qubit purity at equilibrium is given by 

P Eq {e) - 



(7) 



Here the temperature-dependent coupling spectrum 

G? T (w) = G (uj)(n T (uj) + 1) + Go(-w)n T (-w), (8a) 

is written in terms of the average occupation number at 
inverse temperature j3 = 1/T, 



P Eq + e 2 JZ,dwG T (w) 











1 



exp(/3ftw) — 1 ' 
and the zero-temperature bath-coupling spectrum 

Go(w) = e 2 ^2rj 2 k (hw a ) 2 /r 1 2 lax S{Lj - uj k ). 

k 

The equilibrium value purity of the TLS is 
P Eq = tanh(-/3ftw a /2), 



(8b) 



(8c) 



(9) 



with the ground and excited populations respectively 
given by 



Pee =P£ q =(l+P Bq )/2, 
P 99 =P Eq = (l-P E q)/2. 



(10) 
(11) 




FIG. 1: Excitation of the TLS at thermal equilibrium as a 
function of the inverse temperature /3 either with (dashed) 
and without (solid) considering the effect of the system-bath 
interaction. Parameters: memory time of the bath t c = 2/ui a , 
peak of the Lorentzian bath spectrum ojq — 2uo a , maximal 
coupling strength to the bath \rimax\ 2 = w a /100, where hw a 
is the energy separation of the TLS. 

Using a similar analysis, the mean interaction energy 
to 0(e 2 ), is given by 

(HsB^)}Eq = ~hjJa<? 

(13a) 



1 + (2 SZ> dujG T {w) [P+/C+H + P^K-(w)] 
where the quantity in brackets is dimensionless, and 

■ [cosh f3h(uj a =p w) — 1 ± sinh/37i(w a =p w)l . (13b) 
For a Lorentzian coupling spectrum, 



F 2 



?]k — T]max \ 



r 2 + (wo - Lu k y 



(14a) 
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the mean interaction-energy at T = OK, is simply given 
by the bath-induced lamb shift [To) 



{HsB)Eq ~ — fi^a 



duj- 



1 



/„ r 2 + (^ -^) 2 (1 + W/W )' 

(I'll-) 

This proves the negativity of the mean system-bath in- 
teraction energy in equilibrium. 



III. DISTURBANCE OF EQUILIBRIUM BY 
IMPULSIVE QND MEASUREMENT 

The Hamiltonian is intermittently perturbed by the 
coupling of the system (qubit) to the detector (measuring 
apparatus) , designed to effect a QND impulsive measure- 
ment in the <r z -basis. Such a measurement projects the 
qubit onto the |e) or \g) energy states. We stress that the 
measurement results are unread, i.e., the qubit dynamics 
is changed by non-selective measurements. 



measurements), by tracing out the detector degree of 
freedom. The total effect on the system density-operator 



Ps^Tr D {U c ps® \0) DD (0\} 

= \e){e\ps\e)(e\ + \g)(g\p s \g){g\ (20) 

i.e., the diagonal elements are unchanged, and the off- 
diagonals are erased. Since the TLS is entangled with 
the bath, the effect of the measurement in Eq. (fTB)) is: 

Ptot(0) = PEq -+ P f ot = Tr D {UcPtot ® |0) DD (0|} 



= \e)(e\ptot\e)(e\ + \g){g\p t ot\g){g\ 
^ e \e){e\+ P f g \g)(g\. 



(21) 



Since Hsd in Eq. (fH))) commutes with Hs, we 
may consider the measurement-induced evolution of 
(Hsb{t)), rather than (i? tot (r)). In the impulsive limit 
(t — > 0), the measurement yields: 



A. Dynamic description of the measurement 

The time-dependent system-detector coupling (to the 
fcth detector) has the form 

Hsd® = ^(1 + o z ) (|0)(0| + |1)(1| - |0>(1| - |1)(0|) . 

(15) 

where (l + a z ) = |e)(e| ensures QND measurement of the 
qubit energy, and 



h {t) = £. (tanh 2 (LJl | _ [ 



(16) 



is a smooth temporal profile of the system coupling to 
the detector qubits during the measurement that occurs 
at time to and has a duration of t. 

The detector (ancilla) qubits have energy-degenerate 
states |0), |1) so that we may set the detector Hamiltonian 
to be zero 



Hn = 0. 



(17) 



This form of the single- measurement Hamiltonian Hs d 
was chosen so that the measurement interval is [0, r]: 



{Hsb) eq - {Hsb{t)) M 

= Tr {p tot (0) d (0\U&H S b(P)Uc\0) d } • (22) 

Finally, using the RHS of flTjJ} and |@J), we get: 



d(0\uIHsb(0)U c \0) d =0 



H S b(t)) m = 0. (23) 



In fact, this result follows immediately from the nature 
of the projective measurement: 



(H SB ) 



1 



(HsB)Eq 



+ 2 S 7 ^- Tr [( 5 fe + b l) <J ^ a ^PEq} (24) 

k 

= -^{HsB)Eq - -z{H S B)Eq = 0, 

where we have used the identity o~ z o- x o~ z — —a x . 

This expresses the vanishing of Tr {p tot (r)_H sb} M due 
to the diagonality of Ptot( T ) with respect to S. Since 
Hd = 0, the detector mean energy is not affected by the 
measurement. 



e -ifo dtH SD (t)/n _ Uc _ ^ 

where Uc denotes to the CNOT operation. 
In our model (Eqs. (0]) - ([Tg| ) 

= \l) D \e)(e\ + \0) D \g)(g\. (19) 

The measurement consists in letting the TLS interact 
with the detector (a degenerate TLS) via Hg^,. The mea- 
surement outcomes are averaged over (for non-selective 



B. Post-measurement heating 

As shown in (IIII A[) above, a nearly-impulsive (projec- 
tive) quantum measurement (r — > 0) of S, in the \g), \e) 
basis, using the energy supplied by Hsd(0 < t < t) elimi- 
nates the mean system-bath interaction energy. Now the 
pre-measurement equilibrium mean value, {Hsb)eqi IS 
negative, as is shown above (Eq. (|13ap ) by second-order 
perturbation theory, provided the temperature is posi- 
tive, i.e., the \g) state is populated more than the |e) 
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state at thermal equilibrium. Hence 



(Hs B (0)) Eq < ~ (H SB (r)) m = 0, 
(H SD (t)) = -(H SB (t)) M . 



(25) 



After the measurement (as Hso(t > t) = 0), time- 
energy uncertainty at At < l/w a results in the break- 
down of the RWA, i.e., (Hg + H B ) is not conserved as 
At grows. The resulting (Hs) + (H B ) changes stem from 
the non-commutativity of Hsb and H$d- Only (H to t) 
is conserved, by unitarity, until the next measurement. 
Hence, the post-measurement decrease of (Hsb) with At, 
signifying the restoration of equilibrium: 

{Hsb{t)) m = -> {H sb {t + At)) < 0, (26) 

is at the expense of the increase 

(H S + H B ) = (H tot ) - (H SB ) > 0, (27) 

i.e., heating of the system and the bath (Fig. [2J [3]), com- 
bined. 



.x 10" 




FIG. 2: System evolution as a function of time. Excited-level 
population as a function of time for initially zero-temperature 
product state, followed by relaxation to quasi-equilibrium and 
then subjected to a series of measurements (vertical dashed 
lines). Measurements of finite duration (t/j = 0.11/o; a ) (blue 
line) results in somewhat larger heat-up than impulsive mea- 
surements (red line), but the dominant effect is the same 
for both. Observe the agreement between 2 nd order mas- 
ter equation (green), two-quanta exchange with a discrete 
bath, and exact numerical solution for a discrete bath of 40 
modes (black dashed). Parameters: memory time of the bath 
t c = 10/W, peak of the bath spectrum u>o = uj a , maximal 
coupling strength to the bath \rj max \ 2 = 0.07cj a , where hto a is 
the energy separation of the TLS. 



C. Short-time post-measurement qubit evolution 

Let us denote the even part of the bath state by 
|B even )|e) and that of the odd part as \B odd ), then: 




FIG. 3: Excitations as a function of t of the 40 modes in the 
two-quanta model. A is the mode number. Parameters as in 

pig. m 



\B°™(t)) ® Iff) + \K d g(t)) 



!*»,*(*)>• (28) 



Here B even (respectively, B odd ) is a combination of bath 
iV-eigenstates with eigenvalues differing from N by even 
(respectively, odd) numbers. 

The post-measurement evolution of the system alone, 
described by ps = Trsptot, is not at all obvious. Its 
Taylor expansion holds at short evolution times, At 

At 2 

Ps (t + At) ~ Ps (t) + Atpsir) + — p s (r) + . . . (29) 

The 0th order term is unchanged by the measurement, 
Ps(t) = ps(t<0). 

Due to the post-measurement vanishing of the off- 
diagonal elements of p to t (Eq. pij). for ptot(t) — 
|*n, ff (*))<*nrf(*)| (Eq. HMD), wc have 

(Ps) eg (t)=(e\ps(t)\g)=Tr B (e\p to t(t)\g) 

(30 ) 

=(B°™(t)\B° n dd (t)} = 

Hence, ps is diagonal at any time t. 

Its derivative immediately after the measurement, 
P's(t), has the form: 

p s (r) = - ie-^\e)(g\Tr B {B (pf g - pf e )} + H.C. 
=0. (31) 

The same argument goes through upon permuting e «-» 
g everywhere for p tot = |* n)e (i))(* n)e (t)|. 

Hence, the first derivative vanishes at t = r(At = 0) 
due to the definite parity of the bath density-operator 
correlated to \g) or |e). This post-measurement vanish- 
ing of the first derivative, ps(t) = 0, is the condition 
for the quantum Zeno effect (QZE)[l3 0, 0, [H- The 
time evolution of ps is then governed by its second time 
derivative ps(r). 
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For the factorisable thermal state, 

Ptot = z^e-f" 1 * = Z-'e-^-Z^e-^, (32) 
we have: 

pg ={e\p tot \e) = {e\Z^e~^\e)Z^e-^ 
= (ps) e ePB (and e <->- g). 

For this ptot, the second derivative of ps immediately 
after the measurement is (cf. Eq. (|21[) ) 

p s (r) = 2a z Tr B {B 2 (p* g - pg)} . (34) 

The scalar factor is positive: 

= Tt b [b 2 Pb ] ({ps) gg - (ps) ee ) > 0, (35) 

where we have used Tr B {p^ g(ee) } = {ps) gg ( ee ) which 
follows from the definition fEq. pTj) ): Pee{gg) = 
{e(g)\ptot\e(g)). The first factor in (|3"5|) is positive by 
virtue of the positivity of the operator B 2 {B being Her- 
mitian), and the second is positive iff there is no popula- 
tion inversion for the TLS. 

Hence, the second derivative in (|29p is positive shortly 
after the measurement, if there is no initial population 
inversion of the system, i.e., for non-negative tempera- 
ture. 



D. Post-measurement state 



to the ME will be discussed elsewhere. The 2nd order 
ME for ps, on account of its diagonality, can be cast into 
the following population rate equations [llj, dropping the 
subscript S in what follows and setting the measurement 
time to be t = 0: 

Pee(t) = ~ Pgg{t) = R g (t)p gg - R e (t)p ee , (38) 

/oc 
dwG-r (w)sinc [(w =p w a )t] . (39) 
-oo 

Here sinc(a;) = 52i£), \y e shall assume that Gq(uj), the 
zero-temperature coupling spectrum, has peak coupling 
strength at loq and spectral width ~ l/t c . 
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FIG. 4: Relaxation rates, R g (dashed) and R e (solid) as a 
function of time. Parameters as in Fig. [2] 



The combined (system- and bath-) equilibrium state 
satisfies: 

pf ot = Z- X e-^ = Ptot = Z -i e -P{H 0+ o ( H% B ))^ (36) 

Thus, for sufficiently weak coupling, Eq. (|33|) dominates. 

How is this reconciled with the non-unitary nature of 
the projection, whereby the mixedness of the total state 
must increase? Indeed, 

Tr[(/') 2 ] = l -Tv[{p Eq ) 2 ] + \^[{a z p Eq a z p Eq )] (37) 

Yet, in the weak-coupling limit, the increase in mixed- 
ness due to measurement is = 0(e 4 ) and hence can be 
neglected. 

IV. POST-MEASUREMENT FREE EVOLUTION 
OF THE QUBIT 

The evolution of ps at longer times (in the 
regime of weak system-bath coupling) may be approx- 
imately described (as verified by our exact numerical 
simulations ) by the second-order non-Mar kovian mas- 
ter equation (ME) [ijj (Fig. [5]) . Higher-order corrections 
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FIG. 5: R g (t) and R e (t) (Eq. ([39])) depicted as spectral over- 
laps of coupling spectrum (blue solid) and sinc((u ± u) a )t) 
(black dashed). 

The entire dynamics is determined by R e (g) (t) 
(Figs. HI [5]), the relaxation rates of the excited (ground) 
states: 

(i) At short times t <C l/w a *C t c the sine function in 
(|3"9"| is much broader than Gt- The relaxation rates R e 
and R g are then equal at any temperature, indicating the 
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complete breakdown of the RWA discussed above: \g) —> 
\e) and |e) — * \g) transitions do not require quantum 
absorption or emission by the bath, respectively. The 
rates RrUA then become linear in time, manifesting the 
QZeH13,|I|: 

i2 e ( B) (*<tc) ~2R Q t, (40) 

/OO 
duG T (uj) = (B 2 ). (41) 
-oo 

This short-time regime entails the universal Zeno heating 
rate: 

J t (Pee ~ Pgg) ~ 4R t(p gg - Pee ). (42) 

(ii) At intermediate non-Markovian times, t ~ l/u) a , 
when the sine function and Gt in R-def have compara- 
ble widths, the relaxation rates R e ( g ) (t) exhibit several 
unusual phenomena that stem from time-energy uncer- 
tainty. The change in the overlap of the sine and Gt 
functions with time results in damped aperiodic oscilla- 
tions of Re(t) and R g (t), near the frequencies luq — 0J a 
and loq + u) a , respectively. This oscillatory time depen- 
dence that conforms neither to QZE nor to the converse 
AZE of relaxation speedup [3, HE 13, will henceforth 
be dubbed the oscillatory Zeno effect (OZE). Due to the 
negativity of the sine function between its consecutive 
maxima, we can have a negative relaxation rate, which is 
completely forbidden by the RWA. Since sine [(u> + oj a )t] 
is much further shifted from the peak of Gt(oj) than 
sine [(lu — u> a )t], R g (t) is more likely to be negative than 
R e {t) (Figs. |U[5]). Hence, p gg {t) may grow at the expense 
of p ee (t) more than allowed by the thermal-equilibrium 
detailed balance. This may cause transient cooling, as 
detailed below. 

(hi) At long times t 3> t c , the relaxation rates attain 
their Golden- Rule (Markov) values [TlT| 

R e{g) (t »t c ) ~2^G T (±^ a ). (43) 

The populations then approach those of an equilibrium 
Gibbs state whose temperature is equal to that of the 
thermal bath (Fig. EJ). 

If we repeat this procedure often enough, the TLS will 
either increasingly heat up or cool down, upon choosing 
the time intervals Atk to coincide with either peaks or 
troughs of the p ee oscillations, respectively. Since consec- 
utive measurements affect the bath and the system dif- 
ferently, they may acquire different temperatures, which 
then become the initial conditions for subsequent QZE 
heating or OZE cooling, Fig. [6] The results are shown for 
both different and common (Fig. [7|) temperatures of the 
system and the bath. Remarkably, the system may heat 
up solely due to the QZE, although the bath is colder, 
or cool down solely due to the OZE or AZE, although 
the bath is hotter. The bath may undergo changes in 
temperature and entropy too (Fig. [3]) . 
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FIG. 6: Example of a system experiencing first Zeno heating, 
then oscillatory-Zeno cooling, obtained from the second-order 
master equation (black-solid) and from the exact numerical 
solution for a discrete bath of 40 modes (blue dashed). 




FIG. 7: Maximal Zeno heating (blue solid) and subsequent 
maximal cooling (black dashed) as a function of common ini- 
tial temperature of system and bath as — olb = huj a (3- Note 
the critical temperature for oscillatory-Zeno cooling. Param- 
eters: memory time of the bath t c — 10/u> a , peak of the bath 
spectrum uio = w a /0.7, maximal coupling strength to the bath 
\Vmax\ 2 = 4.36a; a . These effects can be strongly magnified by 
choosing other suitable parameters. 



V. DERIVATION OF COOLING CONDITIONS 



By integrating Eq. (|3T)|) over time to acquire J g ^(t), 
one arrives at the following result: 



/•oo 

Jg(t) = / dujGo{u>)nT{to)smc 2 ((u) — u) a )t) 
Jo 

/•oc 

+ / duG (uj)(n T (uj) + l)sinc 2 ((w + uo a )t) (44) 
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Je(t) = / dwGo(w)nr(w)sinc ((u) + u> a )t) 
Jo 



+ dujG Q (uj){n T {u) + l)sincr((a; - ui a )t) (45) 
Jo 

To obtain cooling below the equilibrium temperature, one 
requires that: 



Je(t) > Jf« 



J g (t) J, 



Eq 



R 



Eq 



n{ujg) 
n(uj a ) + 1 



(46) 



Rearranging the terms in the above equation, gives the 
cooling condition, 



GoM- 



sin 



(uj - U> a )] 



> 



G (u>)- 



{U> - LU a ) 2 

2 [l^+^)] 
(uj + UJ a ) 2 



(n T (uj a ) - n T (u)) 

(n T {uJ a ) + n T {uj) + 1) 

(47) 



A general quest for finding the spectral density func- 
tion Go (w), which satisfies the above condition in some 
time interval, at any given temperature T, is quite dif- 
ficult. In the high temperature limit i.e., nr(u>) >> 1 
one can find a necessary condition on the peak position 
of Gq(u>), which can satisfy the above inequality. Sub- 
stituting the high-temperature limit for nr{u>) — l/(3u), 
one can show that in order to allow cooling Gq{ui) needs 
to be concentrated in the frequency interval defined by 



< LU < fl, 



where 



pn 



i 



0w tt + v/4 + + pul 



(48) 



(49) 



Though /30 > 1, it is only the maximum possible bound 
on the detuning of the bath spectrum for the qubit fre- 
quency, indicating that one should not detune the bath 
spectrum too far from uo a to see the cooling effect. Wc 
have numerically verified these conditions for various 
bath coupling spectrums. In the same spirit one can find 
regions in frequency space, where for specific times there 
will be no cooling, independent of the shape of Gq(ui). 



VI. ENTROPY DYNAMICS 

One may always define the entropy of ps relative to its 
equilibrium state po ("entropy distance") and the nega- 
tive of its rate of change, as[3,[^]: 



S(ps(t)\\Po) 



Tr{p S (t) hips(t)} - Tr{p s (t) \np } (50) 



(T {t) = --S(p s (t)\\po). 



In the Markovian realm cr(t) > 0@, [1,0 is a state- 
ment of the second law of thermodynamics. Since 
ps is diagonal, it follows that cr(t) is positive iff 
4r \p ee (t) — (po)ee| 5= 0, consistently with the interpre- 
tation of the relative entropy S(ps\\po) in (|5Tj) as the 
entropic "distance" from equilibrium. Conversely, when- 
ever the oscillatory p ee (t) drifts away from its initial or 
final equilibria, <t takes negative values (Fig. [5]). 
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FIG. 8: a(t) (negative of relative entropy rate of change). 
Parameters as in Fig. [2] 



VII. DISCUSSION: REALIZATION AND 
PRACTICAL CONSEQUENCES 

Consider atoms or molecules in a microwave cavity 



(Fig. 9(a) ) with controllable finite-temperature coupling 
spectrum Gt(^) centered at luq. Measurements can be 
effected on such a TLS ensemble with resonance fre- 
quency uj a in the microwave domain, at time intervals 
Atk ~ 1/(^0 ±w a ), by an optical QND probe jl9( at fre- 
quency uj p (jj a ,ujQ. The probe pulses undergo different 
Kerr-nonlinear phase shifts A0 e or A0 g depending on 
the different symmetries (e.g., angular momenta) of |e) 
and \g). The relative abundance of A</> e and A<t> g would 
then reflect the ratio p ee (tk) / Pgg{tk)- Such QND probing 
may be performed with time-duration much shorter than 
lo^ 1 , i.e. LL> a T~k *C 1, without resolving the energies of |e) 
and \g). 

An experimental scenario involves collective N-atom 
coupling to near-resonant RF resonator (Fig. 9(a) I. Let 
us choose ground sublevels \m = — 1) = \g), \m = +1) 
i . with Zeeman splitting ui eg = MHz. The collective 
Rabi frequency of N > 10 s atoms at a cavity antinode 
is N 1 / 2 Q eg /2'K = lOOKHz. An optical beam will rotate 
in polarization (Fig. |9(b)[ ), thus performing QND mea- 
surement (readout) that resolves \g) and |e) (by their 
symmetry, not by energy) if its Rabi frequency: 



J eg- 



(52) 
(51) 

Such a Rabi frequency corresponds to RWA violation, as 
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(b)Energy level 
description 



(a)Atoms in microwave cavity 

FIG. 9: Possible experimental setup 



discussed in the text above. 

Non-selective measurements increase the Von- 
Neumann entropy of the detector ancillae. Since our 
ancillae are laser pulses, they are only used once and we 
may progressively change the TLS ensemble thermody- 
namics by consecutive pulses, disregarding their entropic 
or energetic price. 



VIII. CONCLUSIONS 

To conclude, we have shown that frequent QND 
measurements may induce either anomalous heating or 
anomalous cooling of TLS coupled to baths on non- 
Markovian time scales. These findings defy the standard 
notions of quantum thermodynamics regarding system 
equilibration in the presence of a thermal bath. 

The practical advantage of the predicted anomalies is 
the possibility of very rapid control of cooling and en- 
tropy, which may be attained after several measurements 
at t > <jJ~ x and is only limited by the measurement rate. 
By contrast, conventional cooling requires much longer 
times, f > i c , to reach thermal equilibrium. 
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